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Density functional theory (DFT) is applied to atomic spectra under perturbations 
of superfluid liquid helium. The atomic DFT of helium is used to obtain the dis- 
tribution of helium atoms around the impurity atom, and the electronic DFT is 
applied to the excitations of the atom, averaging over the ensemble of helium con- 
figurations. The shift and broadening of the D\ and D2 absorption lines are quite 
well reproduced by theory, suggesting that the DFT may be useful for describing 
spectral perturbations in more complex environments. 

1 Spectroscopy in liquid helium 

Spectroscopic measurements of impurity atoms and molecules in^iiperfluid 
helium have been attracting considerable interest in recent years.Era The re- 
pulsive force between an impurity and helium atoms induces a bubble around 
the impurity This leads to a weak perturbation of helium atoms on the 
spectra of impurities. The line shifts and spectral shapes induced by the he- 
lium perturbation provide information on the properties of the bubble in the 
quantum liquid as well as the excited states of the impurity. Since the per- 
turbation is weak, this method also provides a unique tool for spectroscopic 
measurements of atomic clusters at low temperature.0 

Because of its .simplicity, perturbations on alkali-atom lines have been 
studied extensivelyBtm For cesium (Cs) atoms, there are two s-to-p transi- 
tions, the Di (si/2 —> P1/2) and D2 (si/2 ~ * P3/2) lines, both of which are 
blue-shifted and acquire widths in a helium bath. The shifts and widths of 
the two lines are different, and the D2 line has a skewed shape suggesting 
a double-peak structure. These features were first analyzed with a collective 
vibration model of the helium bubble.Ha That model reproduced average peak 
shifts, but gave line widths less than a half of observed ones. A more sophis- 
ticated analysis has been made treating the liquid helium environment by the 
Path-Integral Monte-Carlo method.El However, the method is very costly in 
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computer resources and is difficult to apply to more complex systems. We will 
show that a density functional theory (DFT) together with a statistical treat- 
ment of helium configurations provides a simple and quantitative description 
for the helium perturbations. 



2 Application of DFT to impurity spectroscopy 



2.1 DFT-plus- statistical description of liquid helium 

The energy of liquid helium in the DFT is assumed ta have the form, E 
J drWo( r )> where we adopt the Orsay-Paris functional,' 



WoM = ^ Vv^r) 



l - J dv'p{v)p{v')V u {\v 



1+7 



(1) 

Here, m is the mass of a helium atom and p r is a coarse-grained density, and 
Vlj is a screened Lennard- Jones potential. 

The effect of the impurity was treated by including in Eq. (|l|) a potential 
interaction, Vj(r), between the helium atoms and the impurity, 



W(r) = Wo(r)+Vt(r)p(r). 
We approximate the Vi(r) as a contact interaction, 

t r / \ 27TO 

Vi(r) = p e r, 

rrin 
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where m e is the electron mass and p c (r) is the electron density of the impurity 
which is calculated with the electronic DFT in Sec. 2.2. The scattering length, 
a, is determined from the observed low-energy electron-helium cross section. 

Utilizing the energy functional, E[p] = J drH(r), we calculate the density 
profile of liquid helium, putting the impurity atom at the origin. Minimizing 
the grand potential at zero temperature, fl = E[p(r)] — fiN, leads to a Hartree- 
type equation 



Vi(t) 



(4) 



The equation is solved with the boundary condition that the density go to the 
bulk density po at large r. Results indicate a sharp rise in the helium density 
at r Rs 6 A. This corresponds to the bubble radius. 

We use the p(r) computed above to generate an ensemble of configurations 
of helium atoms as follows. Take a large volume surrounding the alkali atom 
and denote it as V. This volume includes N helium atoms on average, where 
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N is given by J v drp(r) = N. We randomly sample N helium positions in 
V according to the density distribution p(r). This sampling procedure gives 
probability distribution, iu(ri, • • • , rjv) = HfL 1 (p(r i )/N). 

2.2 Helium perturbation on spectra 

Orbital wave functions of valence electrons in impurity, ?/>(r), are calculated 
using DFT with Dirac wave functions and kinetic energy operator. We need 
accurate wave functions at large distances from the atom, which cannot be 
achieved with the traditional LDA functional due to the incorrect orbital 
eigenvalues and the incorrect asymptotic behavior of the potential. These 
problems are diminished with the gradient correction which was designed to 
produce the correct asymptotic behaviour of the potential. 

We use first-order perturbation theory to evaluate the orbital shifts in 
the ensemble of helium configurations r = (r l7 --- , rjv). The same he- 
lium configuration is used for the ground state s\/2 and excited states Pin 
and P3/2, following the Frank-Condon principle. For S1/2 and P1/2 states, 
the energy shifts of the valence electron is then calculated as AE^^r) = 
(2ira/m c ) J^. |-0W(rj)| 2 , where k stands for orbital quantum numbers (£j) 
and either m state may be taken. For p 3 / 2 states, the matrix elements de- 
pend on m and we have to diagonalize a 4 x 4 matrix to get the energy shifts. 
We then obtain two eigenenergies, each of which is doubly degenerate. 

Each helium configuration produces an energy shift and possible splitting 
but the transitions remain sharp. The line broadening comes from the ensem- 
ble average over helium configurations. The line shape of the D\ (si/2 —* P1/2) 
transition is given by 

S Dl {E)= J o!tw{t)8 (e - [ae^'^{t) _ AE<- s ^^{t)^ , (5) 

where E is a shift from the energy position of the free atom. For the D2 
( s i/2 — > P3/2) transition, we have a similar expression but need to add the 
two eigenmodes. 

To calculate line shapes of the Cs D transitions in liquid helium, we evalu- 
ated Eq. (^) by sampling 100 000 helium configurations, generated according 
to the DFT density profiles. The calculated energy shifts are added to the 
observed D lines of free Cs atom (A = 894.9 nm for D\ and 852.7 nm for 
£>2)- Then, the intensity is estimated by counting number of events in bins of 
wavelength AA = 0.1 nm. The obtained intensity spectra are shown in Fig. [l]. 
The D\ line can be well approximated by a single Gaussian, while the D2 line 
has a double-peaked structure. n Xhe calculated line shifts and shapes agree 
with experimental observations.tJel 
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(a) D, excitation spectra 
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Figure 1. (a) Cs D\ excitation spectrum at different helium pressure; P = 9i 10, and 20 
atm. (b) The same as (a) but for D2 excitation spectrum. Experimental data|3 are plotted 
as filled squares. 



3 Conclusion 

We have developed a simple model to describe atomic spectra of impurities 
embedded in the superfhiid helium. Various features in the atomic spectrum 
of Cs, including line shifts, broadening, and skewness, are nicely reproduced 
in our calculation without any adjustable parameters. The model is simple 
enough to apply to more complex chromophores such as molecules and clus- 
ters. Detailed analysis is found in our recent paper .0 
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